Investigation of the magnetic fluctuations in Tb2Sn207 ordered spin ice by high 
resolution energy-resolved neutron scattering. 
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We have studied magnetically frustrated Tb2Sn2 07 by neutron difltraction and high resolution 
energy-resolved neutron scattering. At 0.1 K, we observe short range magnetic correlations with a 
typical scale of 4 A, close to the near neighbor distance between Tb'^"'' ions. This short range order 
coexists with ferromagnetic correlations and long range spin ice order at the scales of 18 and 190 
A, respectively. Spin dynamics was investigated at a time scale down to 10~^s, by energy-resolved 
experiments on a backscattering spectrometer. We observe a freezing of the spin dynamics for all 
length-scales, with a strong slowing down of the spin fluctuations when long range order settles in. 
We discuss the spin fluctuations remaining in the ground state in comparison with previous data 
obtained by muon spectroscopy. 

PACS numbers: 71.27.-(-a, 75.25.-l-z 



Geometrical frustration is expected to favor the on- 
set of unusual types of order such as spin liquid and 
spin ices, showing a large degeneracy of the magnetic 
ground state. In spin ices, it leads to non zero ground 
state entropy akin to that of real icei, peculiar freezing 
dynamics^, and short range magnetic orders^. Appli- 
cation of a magnetic field lifts the ice rule degeneracj«^, 
yielding to excitations which span the entire system^, or 
are akin to magnetic monopoles^, depending on the ori- 
entation of the field with respect to the anisotropy axis. 

Canonical spin ices are observed in pyrochlores 
R2Ti207 (R=Dy or Ho), where the rare earth moments 
occupy a lattice of corner sharing tetrahedra. The sta- 
bilization of the spin ice ground state is related to the 
strong anisotropy of the Dy^^^ or Ho"^"*" ions, whose mo- 
ments are constrained to lie along the <111> local axes 
connecting the center of each tetrahedron to the summits. 

The Tb pyrochlores offer a more complex but even 
richer behavior, due to the smaller anisotropy of the Tb'^+ 
ion, and to the fact that superexchange and dipolar in- 
teractions between near-neighbor Tb'^^ ions nearly com- 
pensate. Tb2Sn207 is an intriguing example of an or- 
dered spin ice*. Contrarily to classical spin ices which 
do not order at large scale, here the four tetrahedra of 
the unit cell are identical, yielding magnetic Bragg peaks 
and long range order, at a length scale which increases 
with decreasing temperature T and reaches about 190 A 
at T=0. The onset of this magnetic order is observed 
at Ti=1.3(l) K. An upturn of the correlation length and 
magnetic moment, together with a peak in the specific 
heat, occurs at Tc=0.87(2)K. 

Magnetic fluctuations play a prominent role in the or- 
dered spin state of Tb2Sn207. They were first evidenced 
by the reduction of the ordered moment as measured us- 
ing the hyperfine Schottky anomaly of the specific heat 
and compared with its value from neutron diffractionS. 
Then they were directly observed using /iSR~il2,. Sur- 
prisingly, no evidence of any static component appears 
in the local field probed by /iSR in spite of the presence 



of the magnetic order. The fluctuations of the local fleld 
occur at a time scale estimated to either to 8 10~^^s or 
to 5 10~^ s, depending mostly on the value assumed for 
the local field^ii^. Under applied field, damped oscilla- 
tions in the muon polarization, and a thermal hysteresis 
in the longitudinal relaxation rate, show the presence of 
local configurations frozen at the time scale of the muon 
probe^". 

Very recently, several inelastic neutron scattering mea- 
surements were performed to study the magnetic fluctua- 
tions in Tb2Sn207. Crystal field excitations^^ were mea- 
sured down to 1.4 K and compared to those in Tb2Ti207, 
which remains spin liquid down to about 50 mK. It was 
shown that the wave functions describing the ground and 
first excited states of the Tb^+ ion are exchanged in the 
two compounds However, these differences in the lo- 
cal states of the Tb'^+ ion cannot explain the different 
cooperative ground states in the Ti and Sn compounds 
within the frame of current theories. 

Neutron spin echo (NSE) measurements^^ were per- 
formed in Tb2Sn207 for a moment transfer Q=0.08 A~^, 
outside the Bragg peak positions, probing the dynamics 
of the ferromagnetic correlations. They clearly observed 
some spin dynamics down to 0.8 K, but failed to detect 
any below. Finally, the magnetic cross section was mea- 
sured by polarized neutrons^'^. The spin dynamics was 
deduced from the comparison between diffraction mea- 
surements (integrating the magnetic signal over all ener- 
gies up to 3 meV) and elastic scattering with a coarse 
energy resolution (integration over an energy window of 
±0.15 meV). At 0.04 K, the authors conclude to the pres- 
ence of static Bragg peaks, together with a liquid-like 
structure factor arising from spins moving faster than 
0.04 THz. 

Here we have investigated the spin correlations and 
dynamics in Tb2Sn2 07 by combining neutron diffrac- 
tion and energy resolved neutron scattering measure- 
ments, with a very high energy resolution of 0.5 fieV at 
half width half maximum, probing time scales up to 1.3 
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10~^ s. Both measurements were performed in the tem- 
perature range 0.1 K-1.4 K. The Q and T dependence 
of the magnetic scattering was analyzed quantitatively, 
considering not only the Bragg scattering, but also the 
medium and short range correlations. Diffraction data 
measured at 0.1 K show the coexistence of several length 
scales, corresponding to first neighbors (^4 A), medium 
(18 A), and long range (190 A) correlations. The spin 
dynamics was investigated in different Q ranges where 
either short or long length scales are dominant. When 
T decreases, we observe a freezing of the spin dynamics 
for all length scales, with a stronger effect in the region 
of the transition. Some spin fluctuations persist down to 
the lowest temperature of 0.1 K. We discuss these results 
with respect to the fluctuations in the ground state in 
Tb2Sn207. 
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FIG. 1: (Colour online) a: Magnetic intensity measure< 
Tb2Sn207 by neutron diffraction at small angles, versus th 
moment transfer Q for several temperatures. A spectrun 
above 1.2 K was subtracted. Solid lines are Lorentzian fit 
(see text). In inset, the total diffraction pattern at 0.1 K 
showing the small angle intensity together with the magneti 
Bragg peaks, b: medium correlation length Lc = 1// 

(full squares) versus temperature T. The solid line is a guid 
to the eye. c: norm A of the Lorentzian (open triangles) an( 
intensity of the magnetic Bragg peaks (full squares) versu 
temperature. 



(ILL) with an incident neutron wavelength of 2.52 A. The 
magnetic long range order was already studied in ref. 
from the same data. The magnetic cross section was iso- 
lated by subtracting a spectrum measured at 1.2 K. The 
magnetic diffraction pattern at 0.1 K (inset Fig la) shows 
an intense small angle neutron scattering (SANS) signal 
for Q <0.2 A-i, arising from ferromagnetic correlations. 
The SANS signal (Fig. la) strongly increases with de- 
creasing T, and flattens in the range of the transition. 
It was fitted by a Lorentzian function /(Q) = 4 k'^+Q'^ • 
Good fits are obtained down to about 0.6 K, but the fit 
quality decreases below. The fit however allows one to 
estimate the medium correlation length Lc '^^^ = 1/k 
(Fig. lb) versus temperature in the ordered region. With 
decreasing temperature, Lc ^'^^ increases from values 
close to the near neighbor distance between Tb^+ mo- 
ments (d=3.686 A), to about 18(6) A below 0.8 K. This 
value is about 10 times smaller than that deduced from 
the width of the magnetic Bragg peaks (Lc = 190 A as 
shown in ref. 8) but its temperature dependence is quite 
similar. The parameter A increases with decreasing tem- 
perature. Its variation compares well with that of the 
LRO squared magnetic moment deduced from the Bragg 
peak intensity, although it is less sharp in the transition 
region (Fig. lc). 
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FIG. 2: (Colour online) a: Magnetic diffraction spectrum of 
Tb2Sn207 versus the moment transfer Q at 0.1 K. A spec- 
trum at 1.2 K was subtracted. The solid line is a refinement 
taking into account both LRO and SRO contributions. A lin- 
ear background is also added. The arrows correspond to the 
Q-bands selected in the time-resolved experiment. White ar- 
rows correspond to the Bragg Q-bands, black arrows to the 
low, medium and high Q-bands.; b; (resp c:) calculated LRO 
(resp. SRO) contribution to the diffraction pattern. 



I. NEUTRON DIFFRACTION 

Neutron diffraction measurements were performed O] 
the DIB spectrometer at the Institut Laiie Langeviu 



The short range order (SRO) yields a diffuse magnetic 
scattering which remains clearly visible at 0.1 K. This 
scattering was not analyzed in ref. d. Here we refined it 
as a magnetic structure, assuming that it has the same 
symmetry as the LRO, but a much shorter correlation 
length. The relative contributions of the SRO and LRO 
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FIG. 3: (Colour online) Inelastic neutron spectra of 
Tb2Sn207 for several temperatures measured in the Q-range 
of the 111 and 200 Bragg peaks. The solid line is a fit as de- 
scribed in the text. Dotted, dashed-dotted and dashed lines 
correspond to the quasielastic, elastic and flat "background" 
terms respectively. 



are shown in Fig. 2. From the fit, we obtain values of 
LRO and SRO moments of 5.8(1) and 3.3(3) /ib respec- 
tively. The value of the LRO moment, ordered at the 
scale of 190 A deduced from the width of the Bragg peaks, 
agree with the value previously determined of 5.9(l)/iB. 
The SRO moment compares well with the value of about 
3 obtained by calibration of the diffuse magnetic scat- 
tering measured with polarized neutrons^'^. This analysis 
shows that a large amount of the Tb^+ moments remains 
ordered at a very short length scale, estimated here to 
4(1) A from the refinement, therefore comparable to the 
distance between first neig hbor Tb3+ moments (3.686 A). 
Assuming that both types of correlations contribute to 
the neutron intensity independently, we estimate the to- 
tal ordered moment to M= 6.6(2) /xb /Tb. We notice 
that this value is exactly the value calculated for the to- 
tal local moment within the crystal field model (Fig. 15 
of ref. [n!). 

The above analysis shows the presence of three char- 
acteristic length scales for the spin correlations at low 
temperature, with values ranging from the first neighbor 
distance to length scales of about 18 cubic unit cells (the 
lattice constant a = 10.426 A). A distribution of length 
scales was previously inferred from an independent anal- 
ysis of the profile peak shape^^. 



II. ENERGY RESOLVED NEUTRON 
SCATTERING 

The analysis above motivates a study of the slow spin 
dynamics at different Q-ranges, where either short or 
long length scales dominate. For this purpose we per- 



formed energy resolved neutron scattering measurements 
on the IN16 backscattering spectrometer of the ILL, with 
an incident neutron wavelength of 6.27 A. This instru- 
ment offers a dynamic range of ± 15 fj,eV for the ob- 
servation of fiuctuations characterized by a quasielastic 
response. Due to the large angular acceptances of the 
backscattering method the wave-vector transfer resolu- 
tion is rather relaxed, anyhow sufficient to discriminate 
regions of interest with or without contributions from 
the Bragg peaks associated with the long-range order- 
ing process. Optionally one can switch off the energy 
transfer analysis to follow only the elastic response within 
the resolution window of 1 /leV at full width half max- 
imum (FWHM). Accordingly one can obtain a temper- 
ature scan of the elastic contribution over the whole Q- 
range and compare intensity in specific Q-bands to the 
temperature variation observed for the Bragg peaks in 
the diffraction experiment. 

The sample was packed in a 1 mm thick flat can in or- 
der to minimize the sample absorption. The transmission 
was evaluated to 0.9. The counting time was 12 hours for 
each temperature. A vanadium sample of the same shape 
was measured to calibrate the angle dependent detection 
efficiency and determine the energy resolution. 

The spectra were corrected for background and absorp- 
tion using spectra of the empty sample holder measured 
in the same conditions, taking into account the specific 
corrections due to the backscattering geometry^'*. Tem- 
perature scan of the elastic resolved scattering was taken 
in the range 0.1< T <1.4 K. The quasielastic spectra 
were recorded in a range 0.1< Q <1.9 A~^. They were 
grouped into five Q-bands, three bands out of the Bragg 
peaks, and two at the Bragg peak positions. The Q- 
bands outside the Bragg peaks are respectively called 
the low Q-band, (Q = 0.55 ± 0.15 A"-'^), the medium 
Q-band (Q = 1.5 ± 0.1 A'^) and the high Q-band 
(Q = 1.85 ± 0.05 A-i). The Q-bands at the Bragg peak 
positions are respectively Qiii+200 = 1.1 ±0.3 A-i and 
Q220 = 1-7 ± 0.1 k-\ the first for the 111 and 200 
together (not well resolved individually) and the second 
for 220 contribution. The 311 Bragg peak is out of the 
accessible Q-range. We have situated these Q-bands in 
the diffraction pattern of the Fig. 2. One sees that the 
Bragg Q-bands mostly probe the long length scales (the 
Bragg intensity dominates), whereas the low, medium 
and high-Q bands probe the shortest lengthscale. We 
also note that the nuclear Bragg contribution on all of 
the above peaks is either zero or negligible, and that the 
overall nuclear incoherent contribution is also small due 
to the small incoherent cross-sections of the elements. By 
fitting the vanadium patterns with a Gaussian function, 
we get an energy resolution at half width half maximum 
(HWHM) of 0.4 /xeV in the Bragg peak and large angle 
regions, and 0.55 iieV in the low angle region. 

Typical energy resolved neutron scattering patterns 
are shown in Fig. 3, in the Q-band covering the 111 and 
200 Bragg peaks. They clearly show a quasi-elastic signal 
in the temperature range 0.8-1.2 K. The quasielastic sig- 
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nal strongly narrows with decreasing temperature, which 
reflects the freezing of the spin dynamics around the 
transition. The neutron intensity was fitted by the sum 
of three contributions. First, an elastic (or resolution- 
limited) signal of intensity Ci, of magnetic origin mostly, 
accounts for the spins which fluctuate at time scales 
longer than the resolution limit. Second, a Lorenzian 
function of intensity C2 and half-width at half maximum 
(HWHM) r accounts for the quasi elastic signal. The 
Lorenzian quasi elastic signal was convoluted by the in- 
strumental resolution. An energy-independent flat level 
" background" Cb accounts for the " fast" fluctuations (or 
short time scales), corresponding to energy widths larger 
than the energy window. Above 1 K, the flat level and the 
rather broad Lorentzian contributions cannot be deter- 
mined independently with good accuracy. So, the value 
of Cb was fixed to the value at 0.9 K. Below 0.6 K, the 
quasielastic linewidth has the same order of magnitude 
as the resolution, and the flat level vanishes. 

The parameters deduced from the fits are plotted ver- 
sus temperature in Fig. 4 in the Q-band of the 111 
and 200 Bragg peaks. Similar results are obtained for 
the Q-band of the 220 peak. The elastic intensity Ci 
increases with decreasing T. Its temperature variation 
scales quite well with that of the Bragg intensity deduced 
from diffraction data. The latter represents a constant 
angle integral of the scattering cross-section over energy 
transfers up to the incident neutron energy (inset Fig. 
4). Just below Tc=0.87(2)K, the intensity measured by 
diffraction increases more sharply with decreasing T than 
the elastic intensity. The differences in the T-dependence 
of the Bragg peak and the elastic resolved intensities can 
be explained by two factors, one being the possibility of 
integrating the "fast" dynamics (represented by the flat 
level in the energy resolved experiment), the other the 
coarser Q-resolution of the time resolved patterns, which 
integrate a contribution from the SRO not included in 
the Bragg peaks of the diffraction patterns. 

The quasielastic intensity C2 shows a broad maximum 
in the transition region, where the contribution from the 
fast ffuctuations starts to enter the energy window. As T 
decreases further, C2 starts to decrease when the spectral 
weight is transferred to the elastic intensity. The energy 
width 2r decreases rapidly down to 0.8 K; a quasielas- 
tic signal of width close to the resolution limit remains 
at 0.09 K. The energy independent flat level Cb also de- 
creases with decreasing temperature. At 0.09 K, there is 
no visible contribution of the "fast" fluctuations to the 
energy pattern and Cb was fixed to zero. 

The temperature variations of the parameters Ci, C2, 
Cb and F reflect a freezing of the spin dynamics. The 
elastic contribution increases at the expense of the " fast" 
fluctuations and quasielastic contributions, whereas the 
decrease of the linewidth reflects the increase of the re- 
laxation time probed in the time window of the measure- 
ment. 

The Lorentzian line shape of the quasielastic signal cor- 
responds to a relaxation process which involves spin fluc- 



tuations with a lifetime r (<S(0).S(t)>cx e"*/'^). Taking 
the Fourier transform of this expression yields t—H/T, 
where ^1=0.658 10"^ /leV. s (see note[ll). The T depen- 
dence of the relaxation time is shown in Fig. 4 c. At 
0.09 K, the time scale of the spin fluctuations is equal to 
To=1.3 10-*^ s (r= 0.5 ^ieV). 

Energy spectra in the low, medium and high Q-bands 
also probe a freezing of the spin dynamics. Here, no 
contribution of a flat level was detected and the fits were 
performed with Ci and C2 contributions only. Their vari- 
ations (Fig. 5) reflect those at the magnetic Bragg peaks. 
The F values at 0.09 K are higher, probing fluctuations 
at shorter time scales (tq ^ 0.3-1 10~^ s). The rela- 
tive weight of the "fluctuating" spins with respect to the 
" static" spins at low temperature is also higher in the Q- 
bands outside the Bragg peaks. Here the terms "static" 
and "fluctuating" refer to the time scale defined by the 
instrumental resolution of IN16 (1.3 10^^ s) and sensi- 
tivity of the fit. 

One notices that due to the limited time window and 
the strong change of the correlations with temperature, 
the opposite variations of the quasielastic and elastic in- 
tensities during the freezing process do not compensate 
in a given Q-band. Namely, at the Bragg peaks the "loss" 
in the quasielastic intensity with decreasing temperature 
is much smaller than the "gain" in the elastic intensity, 
whereas the reverse situation occurs outside the Bragg 
peaks. 



III. DISCUSSION 

The above features show a distribution of time scales, 
which extends to longer times with decreasing tempera- 
ture, as the characteristic length scales of the spin cor- 
relations increase. Somewhat similar effects are seen 
in other frustrated magnets. A distribution of fluctu- 
ation times concomitant with a distribution of correla- 
tion lengths is observed by neutron scattering in disor- 
dered systems such as spin glasse a^^'^^ . A broad range 
of time scales in absence of long-range order but with re- 
markably energy(time)-scale invariant short-ranged cor- 
relations was also seen in the well-known geometrically 
frustrated Kagome spin system SCGO^*. 

The time scale of the fluctuations is in rather good 
agreement with the time scale probed by neutron spin 
echo (NSE) experiments^ in the range of the transition 
(t - 0.2 10~9 s for T=0.8 K and Q=0.08 A"!). At 
0.1 K, when NSE did not observe further dynamics, our 
experiment in the backscattering geometry still suggests 
fluctuations at a time scale tq with tq ranging from 1.3 
to 0.3 10~^s, the shortest time scales corresponding to 
the SRO. Fluctuations at longer time scales than tq do 
exist (the elastic scattering), but no faster fluctuations 
are probed (there is no residual background in the energy 
window) . 

The life time of the spin fluctuations may be compared 
to that probed by /iSR experiments^iiS. In ref. [lO, the 
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/j,SR response was interpreted by a relaxation process 
between different configurations of tfie local field B;oc- 
The typical time Tc of this process can be evaluated from 
the longitudinal relaxation rate Xz by the relation 
Xz / ilfj.^ioc)'^ where 7^ is the muon gyromagnetic ratio. 
At 0.1 K Tc could be evaluated to 8 10"" s in ref. Q (tak- 
ing Bioc^O.2 T) and 5 IQ-^ s in ref. M (taking B;oc=0.02 
T). Our estimation of the life time (tq ^--^0.3-1.3 10~^s) 
situates in between these values. 

On the other hand, the fluctuations at longer time 
scales than tq, clearly seen by neutrons (they give rise to 
the elastic scattering and to a static response in the NSE 
experiment), do not yield any static component in the 
muon spectra. A strong fully static component cannot 
be excluded by the neutron data, but seems to be incom- 
patible with the muon ones. To investigate this point fur- 
ther, we computed the muon depolarization in presence 
of a dipolar field randomly reversing with the frequency 
u, as done in ref. 19. This shows that fluctuations with a 
frequency h'>2j ^Jiioc should wash out the precession sig- 
nal, clear-cut signature of a static field in the muon data. 
For Bzoc=0.02 T, this yields a higher limit of 2 lO'^s 
for the typical time scale t=1/i/ of the fluctuations^^, 
which is still compatible with the lower limit of 1.3 10~^s 
given by the instrumental resolution of IN16. This is no 
longer true for the higher value B;oc=0.2 T. Comparison 
between muon and neutron data therefore favor a lower 
estimation of the internal field. 

The fluctuations of the local field in ordered spin ice 
Tb2Sn207 found by /iSR found by are quite similar to 
those observed in the parent spin hquid Tb2Ti207. In 
both compounds, the short scale fluctuations may involve 
tunnelling or thermally activated excitations between the 
six degenerated conflgurations of the spin ice structure. 
In ordered Tb2Sn207, such excitations may also help to 
"switch" domains of longer length scales. 

In Tb2Sn207, our results suggest a homogeneous pic- 
ture where different correlation lengths involve differ- 



ent spin components and different time scales, rather 
than a coexistence of independent clusters with differ- 
ent sizes. In the pyrochlore spinel ZnCr204 and related 
compounds, a cluster-like scattering process has been in- 
voked, with the presence of emergent excitations within 
loops of strongly bound hexagonal clusters^-. More re- 
cently however, theoretical work on the "true" spin ice 
Dy2Ti2 07 showed that this picture does not necessarily 
imply real independent clusters and that long range in- 
teractions are relevant^!. Whether such picture can be 
extended to describe excitations in the ground state of 
the ordered spin ice Tb2Sn2 07, in which we observed 
multiple length scales, remains a matter for further the- 
ories. 



IV. CONCLUSION 

In conclusion, in Tb2Sn207 ordered spin ice we have 
observed a freezing of the spin fluctuations which involves 
the coexistence of different correlations lengths and times 
scales. The distribution extends to longer times with 
decreasing temperature. A marked slowing down is ob- 
served in the transition region, in agreement with pre- 
vious muon and neutron spin echo experiments. At the 
lowest temperature, the time scale of the spin fluctua- 
tions ranges between 1.3 and 0.3 10~^ s. Fluctuations at 
longer time scales are also observed by neutrons, whereas 
they give no clear-cut signature in the muon spectra. 
This may be compatible with the low estimation of the 
local field (B/oc=0.02 T) found in /iSR. The presence of 
multiple time scales remains to be understood in further 
theoretical work. 

We thank B. Frick and R. Chung for the help with the 
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help with the neutron diffraction experiment, and J.-L. 
Ragazzoni for setting up the dilution equipment. 



^ A. P. Ramirez, A. Hayashi, R.J. Cava, R. Siddhartan, B. 

S. Shastry, Nature 399, 333 (1999). 
^ J. Snyder, J.S. Slusky, R.J. Cava, R Schiffer, Nature 413, 

48 (2001). 

3 S.T. Bramwell, M.J. Harris, B.C. den Hertog, M.J.R Gin- 
gras, J.S. Cardner, D.F. McMorrow, A.R. Wildes, A. Cor- 
nelius, J.D.M. Champion, R.G. Melko, T. Fennell Pliys. 
Rev. Lett. 87, 047205 (2001). 

* S.T. Bramwell, M.J. P. Gingras Science 294, 1495 (2001). 
^ H. Fukazawa, R.G. Melko, R. Higashinaka, Y. Maeno, 

M.J. P. Gingras, Phys. Rev. B. 65, 54410 (2002). 
" L.D.C. Jaubert, J.T. Chalker, P.C.W. Holdsworth, R. 
Moessner Phys. Rev. Lett. 100, 67207 (2008). 
C. Castelnovo, R. Moessner, S.L. Sondhi, Nature 451, 42 
(2008). 

* L Mirebeau, A. Apetrei, J. Rodriguez- Car vajal, P. 
Bonville, A. Forget, D. Colson, V. Glazkov, J. P. Sanchez, 
O. Isnard and E. Suard, Phys. Rev. Lett. 94, 246402 



(2005). 

^ P. Dalmas de Reotier, A. Yaouanc, L. Keller, A. Cervel- 
lino, B. RoessU, C. Baines, A. Forget, C. Vaju, P.C.M. 
Gubbens, A. Amato, P.J.C. King, Phys. Rev. Lett. 96, 
127202 (2006). 

° F. Bert, P. Mendels, A. Olariu, N. Blanchard, G. CoUin, 
A. Amato, C. Baines and A.D. Hillier, Phys. Rev. Lett. 
97, 117203 (2006). 

^ I. Mirebeau, P. Bonville, M. Hennion Phys. Rev. B. 76, 
184436 (2007). 

^ Y. Chapuis, A. Yaouanc, P. Dalmas de Reotier, S. Pouget, 

P. Fouquet, A. Cervellino, A. Forget, J. Phys. Cond. Mat. 

19, 446206 (2007). 
^ K.C Rule, G. Ehlers, J.R. Stewart, A.L. Cornelius, P.P. 

Been, Y. Qiu, C. R. Wiebe, J.A. Janik, H.D. Zhou, D. 

Antonio, B.W. Woytko, J. P. Ruff, H.A. Dabkowska, B.D. 

Gaulin and J.S. Gardner Phys. Rev. B. 76, 212405 (2007). 
* O.G. Randl, ILL report 96RA07T (1996). 



6 



The equation which yields the typical time scale holds be- 
cause it involves an energy measurement and a Fourier 
transform. This is the case for Quasielastic or Inelastic 
Neutron Scattering, Electronic Paramagnetic Resonance, 
Mossbauer spectroscopy or energy-resolved Nuclear Mag- 
netic Resonance (NMR). The time involved results from 
a dynamical broadening or from the peak lineshape as- 
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l/tj=l/27r!/, where u is the frequency equivalent to the en- 
ergy fiw involved in the process. On the reverse, when the 
typical time scale is determined from a time measurement 
such as muon spectroscopy or time-resolved NMR, where 
the process involved is expressed as e~*^'^, then t—1/v is 
expressed as the inverse of the frequency. This is also the 
case when r is simply defined with respect to the frequency 



of a physical event, as being the average time between two 
events. 
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FIG. 4: (Colour online) Temperature dependence of the pa- 
rameters deduced from the fit of the inelastic spectra in the 
region of the 111 and 200 Bragg peaks, a: elastic intensity 
Ci and quasielastic intensity C2 deduced from the fit of the 
energy patterns. In inset: the elastic intensity (full dots) is 
compared to the intensity from eleistic scans (open squares) 
and to the integrated intensity of a magnetic Bragg peak mea- 
sured in diffraction (diamonds). All quantities are scaled to 
vary between 1 at 0.1 K and at 1.4 K. b: quasielastic width 
2r versus temperature. In inset: the background level Cb, 
involving fluctuations at a shorter time scale than probed in 
the energy window, c: the fluctuation time T=h/T. Solid 
lines are guides to the eye. 




FIG. 5: (Colour online) Temperature dependence of the in- 
tegrated elastic and quasi elastic intensities. The background 
level was fixed to zero at all temperatures. Solid lines are 
guides to the eye. Fig. 5a, 5b, and 5c correspond to low, 
medium and high Q-bands respectively. In the insets : the 
quasielastic width 2r versus temperature. The width at 0.09 
K in the inset of Fig. 5a was fixed. 



